Causal association of electromagnetic signals using the Cayley Menger determinant 
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In complex electromagnetic environments it can often be difficult to determine whether signals 
received by an antenna array emanated from the same source. The failure to appropriately assign 
signal reception events to the correct emission event makes accurate localization of the signal source 
impossible. In this paper we show that as the received signal events must lie on the light-cone 
of the emission event the CayleyMenger determinate calculated from using the light-cone geodesic 
distances between received signals must be zero. This result enables us to construct an algorithm 
for sorting received signals into groups corresponding to the same far-field emission. 

PACS numbers: 89.20.Dd, 89.70.-a, 89.90.+n, 02.40.Ky, 02.40.Dr 



The deinterleaving of radar pulses is vital for the suc- 
cessful operation of radar warning receivers [ij. Modern 
radars can emit up to a million pulses per second, hence 
in a multi-emitter environment deinterleaving becomes a 
significant problem. Deinterleaving can be simplified if 
radars are distinguishable on the basis of other known 
pulse parameters such as transmitted power, carrier fre- 
quency or pulse width, however in many situations the 
radar pulse parameters are unknown or vary in an unpre- 
dictable way. There are a number of algorithms that are 
quite successful at deinterleaving pulse trains if the pulse 
repetition interval (PRI) is constant and a sufficiently 
long set of adjacent pulses is recorded 0, Mod- 
ern civilian marine radar however often have a pseudo- 
random PRI pattern to minimise the possible effects of 
interference; hence the techniques described in the above 
references will be ineffectual in deinterleaving the re- 
ceived radar pulses. In this paper we propose an alterna- 
tive method for deinterleaving received radar pulses that 
requires no knowledge of the emitter radar characteristics 
and can be applied to only several pulses; this method is 
based on the concept of the light-cone used to describe 
events in Minkowski space-time . 

While we use radar pulses as an example in this letter 
it is important to realise that the proposed algorithm is 
valid for associating any events that lie on the same light- 
cone. Furthermore although the following treatment as- 
sumes four receivers in 2-f 1 space-time it can be extended 
trivially to five receivers in 3 -I- 1 space-time. 



The Cayley-Menger determinant for four points in Eu- 
clidean space is defined as 



D = det 
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(1) 



where sfj is the squared Euclidean distance between any 
two points in the space. 

It has previously been noted that the Cayley-Menger 
determinant can be used to determine if a surface is 
flat 8] . The Cayley-Menger determinate for a flat surface 
is zero if the distance between points on the surface is the 
geodesic distance which, in general, is not the Euclidean 
distance. The proof of this more general result can be 
verified by following the proof in Euclidean geometry as 
presented in section 3.1 of and realising that the result 
generalises so long as there exists a coordinate system in 
which the metric tensor has constant coefficients, which 
is an equivalent definition of the flatness of a surface . 
It is immediately apparent from this definition that the 
Cayley-Menger determinant is zero in 2 -I- 1 Minkowski 
space-time if the space-time events lie on a plane. If the 
Minkowski metric is used to calculate the space-time in- 
terval then sfj may be positive or negative depending on 
the signature of the Minkowski metric. 

As the Cayley-Menger determinant specified by ([1]) 
only contains terms that are cubic in sf^, for errorless 
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measurements 



points in coordinates defined by ^ is 



^ I If all points lie on a flat surface 

1 O {sfj) If any points do not lie on a flat surface 

(2) 

In 2 + 1 space-time the light-cone is a two dimensional 
surface defined by all the possible paths of a photon emit- 
ted at a particular space-time point. Using the condition 
stated in ([2|) we conclude that the Cayley-Mcnger deter- 
minant is zero for all space-time events that lie on the 
same light-cone and hence are causally connected to the 
same emission event. 

Evaluation of the Cayley-Menger determinant us- 
ing IT]) requires knowledge of the geodesic distance be- 
tween points on the light-cone. Calculation of the 
geodesic distance requires determining the metric ten- 
sor for the embedded conical surface and integrating the 
corresponding geodesic equation [6i] . This process can be 
simplified considerably by choosing an appropriate set of 
coordinates {^,ip,C} so that the metric tensor is diago- 
nal and constant on the surface of the light-cone. If this 
can be done then the geodesic distance is the Euclidean 
distance in these coordinates. 

One such set of coordinates is given by 



^ — r cos ( 
— r sin ( 
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and the corresponding inverse relations 



^arctanf^ 
smC 



= c 



(3a) 
(3b) 
(3c) 



(4a) 
(4b) 
(4c) 



where {r, 9, 0} are the spherical polar coordinates and 
are related to the 2 -t- 1 cartesian space-time coordinates 

by 



r = y/{x - xo)^ + iy- yaV + - ^o)^ (5a) 
'y-y^'^ (5b) 

(5c) 



arctan 



9 = arctan 



X — Xo 



t~tn 



^/{x~ xoY + {y- yoY 



where time is measured in natural units, i.e., c — 1 and 
[to,xo^yo\ is the space-time coordinate of the emission 
event. Note that the coordinates defined by ([3|) are not 
the conical coordinates as defined by [l^ , they have been 
formulated so that the metric tensor on the cone is diag- 
onal and constant. 

The infinitesimal Euclidean distance between two 



ds^ = di^ + di^^ + (C^ -I- + arctan cot j d(^ 

-1-2-0 arctan ^"^^ cot C,d^dC, 

—2^ arctan ^'^^ cot (dipdC 

If we constrain the path between any two points to be on 
the light-cone so that C = f then the geodesic distance 
between any two points = f } {^j, ipj, Q = 

f}is 



(6) 



The ultimate aim of many signal association algo- 
rithms is to provide information so that the emitter may 
be localised. Calculation of the Cayley-Menger determi- 
nant ^ using ^ requires the space-time coordinates of 
the emission event for insertion into ^ and hence would 
not appear to be useful. 

In the far-field limit the geodesic distance on the light- 
cone can be approximated by the Euclidean distance i.e., 

- ^jf+ii'i - i^jf ~ {xi - xjf+{yi - yjf+{U - tjf 

The far field approximation is valid if the difference be- 
tween the geodesic distance on the light-cone and the 
Euclidean distance is much less than the Euclidean dis- 
tance. It can be shown that this condition is equivalent to 
the inter-antenna spatial distances being much less than 
the spatial distance to the emitter, i.e.. 



V [xt - xoY + {yi - yaf > ^ {x^ - x^f- + (y, - y^Y 

as in a typical radar application. 

In realistic scenarios the measured signal time of ar- 
rival is noisy which means that the expected value of D 
from ([T]) will not be zero, even if all the points lie on a 
fiat surface. To calculate the expected value of D in the 
presence of noise we would add a noise term e* to each 
of the receive events, calculate the square distance as a 
function of the noiseless distances and the noise, substi- 
tute these distances into ([l} and obtain an expression for 
the Cayley-Menger determinant in terms of the noiseless 
space-time distances and et. 

This rather complex procedure can be avoided and the 
effect of noise on ([T]) can be approximated by a rather 
simple expression once it is realised the Cayley-Menger 
determinant contains only terms of the type sfjS^;S^„ 
and hence, to lowest order, the effect of noise can be ap- 
proximated by s^CTj , where is the maximum square 
interval between events and cr^ is the variance in the time 
of arrival noise. Using this result it is possible to con- 
struct a hypothesis test: 

i/o : AH recieve events are caused by the same 
emission event 
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The decision rule we use to test between this hypothesis 
is 



{true 
false 



If 
If 



D < 2(imatfsl 
D > 2(1.96crt 



(7) 



The factors 2 and 1.96 come from that fact that the time 
difference variance is twice the time variance and that 
95% of time of arrival measurements lie between ±1.96(Tt 
of the noiseless time of arrival. 

An algorithm for associating space-time measurements 
to a particular far-field emission event in 2 -|- 1 space-time 
is as follows: 

1. Select four signal reception events from a list of 
possibilities 

2. Determine the Euclidean space-time interval be- 
tween all combinations of event quartets 

3. Calculate the Cayley-Menger determinant us- 
ing ([1]). Use this to test Hq using ([7]). If Hq is true 
then the events chosen in step [1] are associated to 
the same emission event. 

4. Repeat steps [1] to [3] until combinations are ex- 
hausted 

As an example of how this algorithm would work in 
practice consider a four element antenna array with four 
antennas placed at the four corners of a three metre 
square, i.e., 



two emissions events whose corresponding receive events 
X = [ct, X, y] are 



xi = [4.0133, 3, of 

X2 = [3.2334, 0,3f 

a;3 = [7.3851, -3, of 

a;4 = [8.1967, 0,-3f 



and 



xi = [7.9263, 3, of 

X2 = [5.4521, 0,3f 

a;3 = [2.0474, -3, of 

Xi = [4.4990, 0,-3f 



Each quartet of events can used to estimate the location 
emitter using the time difference of arrival technique [ll[ . 

The Defence Science and Technology Organisation 
(DSTO) of Australia is building an antenna array sys- 
tem to test the association algorithm outlined in this 
paper. The antenna array will be constructed using a 
precise time interval measuring unit such as the ATMD- 
GPX [l3| ■ Once the evaluation system is built trials will 
be conducted with two and three radars placed at about 
5km from the antenna array. The tests will be conducted 
at the DSTO signals testing facility in Adelaide, Aus- 
tralia. Results of the trial will be reported elsewhere. 
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(8) 



Consider a situation in which three signals arrive at each 
of the four antennas so that 



h = {3.2400,4.0133,7.9263} (9a) 

t2 = {3.2334,3.7789,5.4521} (9b) 

h = {2.0474,6.3819,7.3851} (9c) 

ti = {4.4990,5.3704,8.1967} (9d) 

How many emissions are there and what times are asso- 
ciated with each emission? 

The times given in ([9]) were generated by adding zero 
mean Gaussian noise with 50 picosecond standard de- 
viation to the arrival time of signals generated by two 
emitters located 5km from the centre of the antenna ar- 
ray. The third set of arrival times was generated by a 
random number generator and could represent spurious 
measurements. There are S'* ways in which the sam- 
ple data can be grouped into quartets, however of these 
81 possible combinations only two satisfy the hypothe- 
sis test as stated in ([7]). Using the algorithm outhned 
above we are able to correctly conclude that there are 
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